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ABSTRACT 

We present polarization observations of the 6.7-GHz methanol masers around the massive 
protostar Cepheus A HW2 and its associated disc. The data were taken with the Multi-Element 
Radio Linked Interferometer Network. The maser polarization is used to determine the full 
three-dimensional magnetic field structure around Cepheus A HW2. The observations suggest 
that the masers probe the large scale magnetic field and not isolated pockets of a compressed 
field. We find that the magnetic field is predominantly aligned along the protostellar outflow 
and perpendicular to the molecular and dust disc. From the three-dimensional magnetic field 
orientation and measurements of the magnetic field strength along the line of sight, we are 
able to determine that the high density material, in which the masers occurs, is threaded by 
a large scale magnetic field of ~ 23 mG. This indicates that the protostellar environment 
at ~ 1000 AU from Cepheus A HW2 is slightly supercritical {A ss 1.7) and the relation 
between density and magnetic field is consistent with collapse along the magnetic field lines. 
Thus, the observations indicate that the magnetic field likely regulates accretion onto the disc. 
The magnetic field dominates the turbulent energies by approximately a factor of three and 
is sufficiently strong to be the crucial component stabilizing the massive accretion disc and 
sustaining the high accretion rates needed during massive star-formation. 

Key words: stars: formation - masers: methanol - polarization - magnetic fields - ISM: indi- 
vidual: Cepheus A 



1 INTRODUCTION petitive accretion theo ry JBonnell et alj|l998l) . In the third sce- 

nario, core accretion (Mc Kee & TarJ 120030 . massive stars form 



Massive stars are short lived, yet they dominate galaxy evolution 
due to their strong radiation while enriching the interstellar medium 



through gravitational collapse, which involves disc-assisted accre- 
tion to overcome radiation pressure. This scen ario is similar t o 

with heavy elements when they explode in a supernova. Massive - , . . c . ~ nr. —ni.^j. 

. . the favored picture of low-mass star-formation (Sh u et all 1 19951) . 



stars typically form in distant, dense clusters. In such regions, grav 
itational, radiation and turbulent energies are different from those 



in which magnetic fields are thought to play an important role by 

removing excess angular momentum, thereby allowing accretion 

in the regions that form lower mass stars such as the Sun. Specif- , rr; ~r\ |,„„ J rr „ ' : n 

to c ontinue onto the star (Shuetal. 1995; Basu & Mouschovias 
ically, young stars more massive than 8 solar masses (M a ) have a | , nr . x TT , . . . „ 

, 1994). However, the accretion rates during massive star-tormation 

radiation pressure that would be sufficient to halt infall and prevent ... tl , . , ^ , . r .. „ , 

r I — : ij -i are significantly higher than during the formation of a low-mass 

the accretion of additional mass (e.g. Wolfire & Cassinelh 1987). , ■ ■ , , , ■ ,■ , , 

star, and it is unclear whether massive accretion discs are able 

Still, stars more massive than 200 M n have been observed (e.g. . , . , , t . „ , , 

i— : II 1 . to sustain these rates without support by strong magnetic fields 

| FigeretalJ | 1998|) , indicating that the processes governing massive ( Mc Kee & ^j l^ Stm> if such magnetized dlscs exist around 

star-formation are yet poorly understood. . a ... , 

%. „ massive pr otostars, outflows will arise as a natural consequence 

A number of different formation scenarios have been pro- rr . ' — , • j-„„4 

r— ; — I -i r (Banenee & Pudntz 2007). 

posed (Zinnecker & Yorke 2007). These include formation through 

the merger of less massive stars jBonnell & B"ate]|2005h or through Current observations of magnetic fields in massive star for- 

me accretion of unbound gas from the molecular cloud, the com- mation regions are often limited to low density regions (hydro- 
gen number density rc H2 < 10 6 cnr 3 ) and/or envelopes at scales 
of several 1000 AU. Linear polarization observations of dust also 
* E-mail: wouter@astro.uni-bonn.de only provide information on the magnetic field in the plane of 
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the sky, so they have been yet unable to probe the strength and 
the full structure of the magnetic field close to the protostar and 
around protostellar discs. Recent high-angular-resolution submil- 
limeter dust polarization observations of the massive hot molec- 
ular core G31. 41+0.31 indicate that the gravitational collapse 
shows signs of being control led by the magnetic field scales of 
~ 5000 AU dGirart et al]2009l) . Currently, the only probes of mag- 
netic fields in the high density regions close to the massive pro- 
tostars are masers. Because of their compactness and high bright- 
ness, masers are specifically suited for high angular resolution ob- 
servations of weak linear and circular polarization and can thus 
be used to determine th e magnetic field str ength and morphology 
down to AU scales (see IVlemmingsll200l for a review of maser 
polarization ). Most observations have focused on OH and H2O 
masers (e.g. lBartkiewicz et alj|2005l ; I Vlemmings et al.i r2006a). but 
recently, methanol masers h ave been shown to be excellent probes 
of th e magnetic field (e .g. Vlemmin gs et al.l l2006bl : I Vlem mings 
120081 : ISurcis et al.l l2009h - A number of both OH and methanol 
maser observations indicate that, although individual maser fea- 
tures are only a few AU in size, they pr obe the large scale magnetic 
field around the massive protostar (e.g. lSurcis et alj|2009h . 



We observed Cep heus A, which, at a distance of 700 pc 
dMoscadeUi et al.l |2009|) . is one of the closest regions of active 
massive star-formation. Located in the Cepheus OB3 complex, it 
hosts a powerful extended bipolar molecular outflow that likely 
originates f rom HW2, the brightest r adio continuum source in 
the region faughes & WouterlooJll984h . HW2 is thought to be 
a young protostar of spect ral type B0. 5 with a mass of ~20 M Q 
( Jimenez-S erra et al]|2007l) . The extended outflow appears to be 
driven by a small scale (~700 AU) ther mal radio jet, with ionized 
gas ejected at a velocity of ~500 km s' dCuriel et alj|200^) . In ad- 
dition to the bipolar outflow, Cepheus A HW2 is surrounded by a 
rotatin g disc of dust and molecular gas oriented perpendicular to 
the je t dJimenez-Serra et alj|2007l : iPatel et al]|2005l ; iTorrelles et all 
1 20071) . The jet and disc structure supports the picture where mas- 
sive stars form through disc accretion, similar to low-mass stars. 
The environment of Cepheus A HW2 however, is significantly 
more complex than that of typical low-mass protostars. Besides 
HW2, at least three additional young stellar objects are located 
within an area of ~500 x 500 AU. All thes e sources seem to be 
located within or close to the molecular disc 1 Torrefies et al.ll200~7l ; 
IComito et alj 120071 : Q tmenez- SerraetalJ f2009). The presence of 
OH, H 2 and CH3OH (methanol) masers in or near to the disc of 
Cepheus A HW2 gives us an unique opportunity to study its mag- 
netic field ( Vlemmin gs et al . 2006a; Bartkiewicz et al ] l200l) . with 
recent methanol maser Zeeman-splitting measurements revealing a 
dynamically important magn etic field with a ~ 8 mG component 
along the maser line-of-sight dVlem mings 200j|. 



In §|2 we present the methanol maser observations and in §[3] 
we shortly discuss the radiative transfer tool used to infer the full 
3-dimensional magnetic field geometry. §[4] presents the results of 
our observations and in §[5] we discuss the derived 3-dimensional 
magnetic field morphology and what its implication for the role of 
the magnetic field during high-mass star-formation. Additionally, 
we highlight the potential for methanol maser polarization observa- 
tions of the upgraded Multi-Element Radio Linked Interferometer 
Network (e-MERLIN). 
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Figure 1. Total intensity (I) map of the 6.7-GHz methanol masers around 
Cepheus A HW2. Contours are drawn at 5, 10, 20, 40 & 80 Jy beam" 1 . The 
labels identify the features presented in Tables[T]and|2] The vectors indicate 
the observed polarization angle x- Also indicated is the best fit ellipse to the 
methanol maser feature distribution observed with MERLIN. 



2 OBSERVATIONS 

The 5i - 60 A + maser transition of methanol at 6.7-GHz was ob- 
served on 2006 December 2-4 using 6 of the MERLIN telescopes. 
Total on source time was ~ 27 hr, which was interleaved with ob- 
servations of the phase-calibrator source 2300+638. The resulting 
beam-size was 40x30 mas. We used a 250 kHz bandwidth with 256 
channels for a total velocity coverage of ~ 1 1 km s~' centered on 
the source velocity Vlsr = -3 km s _1 . This provided a velocity res- 
olution of 0.044 km s -1 . To obtain the linear polarization, the data 
were correlated with the full Stokes parameters. For calibration pur- 
poses, the continuum calibrators were observed with the 16 MHz 
wide-band mode. Both 3C84 and 3C286 were also observed in the 
narrow-band spectral line configuration and were used to determine 
the phase off-set between the wide- and narrow-band setup. Flux 
and bandpass calibration was done using the primary flux calibra- 
tor 3C286. Instrumental feed polarization was determined using the 
unpolarized calibrator 3C84 and the polarization angle was cali- 
brated using 3C286. The residual instrumental linear polarization 
was determined to be less than < 0. 1 percent from the observations 
of 3C84 and 2300+638. Self-calibration, in RCP and LCP sepa- 
rately, was performed on the strongest isolated maser feature. After 
calibration, the antenna contributions were re-weighted according 
to their sensitivity at 5 GHz and their individual efficiency. Finally, 
5.12 x 5.12" image cubes were created in Stokes I, Q, and U. The 
rms noise in the emission free channels was ~ 6 mjy beam -1 for 
Stokes I and ~ 9 mjy beam -1 for Stokes Q and U. For the channels 
with strong maser emission we are limited to a dynamic range of 
~ 850 due to residual calibration errors. Because of the dynamic 
range limit we were unable to determine the circular polarization 
due to Zeeman splitting in these observations. The non-detection of 
the ~ 0.17 percent circular polarization observed for the 6.7-GHz 
metha nol masers using the Effelsberg 100-m telescope Iviemmingsl 
2008) is unsurprising considering the rms noise limits of our circu- 
lar polarization observations (lcr ~ 0.12 percent). 



3 MASER POLARIZATION RADIATIVE TRANSFER 

Maser theory describes that the fractional linear polarization P\ of 
the non-paramagnetic maser species such as H 2 0, SiO and CH 3 OH, 
depends on the degree of maser saturation and the angle 6 between 
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the maser propagation direction and the magnetic field. Addition- 
ally, the relation between the measured polarization angle x an d the 
magnetic field angle on the sky <f> B also depends on 9, with the po- 
larization vectors generally parallel to the field for 9 < 9 cr j t ~ 55 ° 
and perpendicular to the field for 9 > 9 a i t (Goldr eichetaill 19731) . 
However, non-magnetic effects, such as full maser saturation and 
anisotropic maser pumping, can generate enhanced linear polar- 
ization and ca n rotate the polarization v ectors with respect to the 
magnetic field (Deguchi & Watson 1990). Thus, before we can use 
polarization observations to determine the magnetic field morphol- 
ogy, we need to estimate the level at which these effects influence 
our measurements. The rotation and generation of linear polar- 
ization depends on the amount of Zeeman- splitting and the level 
of maser saturation described by the Zeeman frequency shift gCl, 
the maser decay rate T and the stimulated emission rate R. When 
gQ. » r > R, the linear polarization direction is determined by the 
magnetic field, even in the presence of anisotropic pumping. When 
the effect of anisotropi c pumping is n ot significant, even gCl > T is 
a sufficient condition dWatsonlllool) . The maser stimulated emis- 
sion rate is given by: 

R = AkT b AQ,/47rhv. (1) 

Here A is the Einstein coefficient for the maser transition, which 
is equal to 0.1532 x 10~ 8 s -1 , and k and h are the Boltzmann and 
Planck constants respectively. The maser frequency is denoted by v, 
and T b and AO. are the maser brightness temperature and beaming 
solid angle respectively. We measure T b < 10 10 K for our masers 
assuming that they are almost resolved with the MERLIN beam of 
~ 30 mas. However, when observed at even higher angular resolu- 
tion, more typical sizes are of the order of a few mas and b right- 
ness temperatures reach of order 10 12 K dMenten et al.ll 1992b . The 
beaming angle AQ. cannot easily be measured directly. Still, those 
of H? Q masers in star-forming reg ions are shown to be AQ. ~ 10" 4 - 
10~ 5 dNedoluha & Watsonl Il99ll) . while those in the envelopes 
of ev olved stars are of order 10~ 2 jVlemmings & van Langeveldd 
2005). It is likely that methanol masers have similar beaming an- 
gles. Taking Af2 ~ 10~ 2 as a conservative upper limit, the maser 
stimulated emission rate R < 0.04 s _1 . 

The maser decay rate T is more difficult to determine. 
For maser transitions in volved infrared pumping, whi ch includes 
methanol, T ~ 1 s -1 dScappaticci & Watsonl 1l992h . However, 
infrared trapping can reduce this value, although likely n ot by 
more than an order of magnitude dGoldreich & Kwanlll974l) . Al- 
ternatively, when the maser transition involves collisions, T = 
0.1n# 2 /10 9 s -1 , with n Hn the hydrogen number density in cirr 3 . 
For the methanol masers this would imply T < 0.1 s -1 . However, 
as the collisional decay rate is smaller than the radiative decay rate, 
the radiative decay rate determined T. We thus assume T = 1 s -1 , 
thoug h other autho rs have assumed the lower value of T = 0.6 s" 1 
dMinier et ai]|2002h . 

We have thus shown that even with conservative estimates for 
R and T, we are in the regime where F > R. However, we need to 
asses if gQ. is sufficiently much larger than R and T so that, even 
in the case of anisotropic pumping, the polarization vectors can 
be used to determine the m agnetic field direction. The analysis by 
iNedoluha & Watsonl ( fl99bT) indicates that the effect of anisotropic 
pumping decreases significantly for transitions with high angular 
momentum, such as the 6.7-GHz methanol masers. This is fur- 
ther confirmed by the fact that we do not observe any relation be- 
tween the maser total intensity and fractional polarization, which 
would be expected when non-magnetic effects contribute to the po- 
larization. However, even without having to consider anisotropic 



pumping, the quantization axis of the maser is potentially not suf- 
ficiently determined by the magnetic field when gQ is not sig- 
nificantly larger than R and T. For a magnetic field strength of 
> 8 mG in the Cepheus A methanol maser region ( Iviemrningsl 
2008), gQ > 10 s" 1 . Although the uncertainty on the magnetic field 
determination was estimated to be of order 25 percent, the limit of 
8 mG assumes a magnetic field along the maser line of sight and 
the tr ue magnetic field, and cons equently gSl, will thus likely be 
larger. lNedoluha & Watsorj dl990l) found that, for a ratio R/T < 0. 1 
and gQ/F > 10, the relation between polarization angle and mag- 
netic field direction hold to within ~ 15°. This improves quickly 
to better than ~ 5° when gQ/T becomes larger. We are thus in the 
regime where the polarization vectors directly determine the mag- 
netic field direction, but adopt an additional 5° uncertainty in our 
determinations of the magnetic field direction on the plane of the 
sky. 

For HtO masers it has been shown to be possible to esti- 
mate the saturation level, as determined by the emerging bright- 
ness temperature T b AQ, and the intrinsic maser line-width AV,- by 
fitti ng both the total intensit y and linear polarization maser spec- 
tra dVlemmings et alj|2006al) . As T b AQ determines the relation be- 
tween Pi and 9 we are then able to constrain 9 using the measured 
Pi and thus provide the full 3-dimensional magnetic field direction. 
We have adapted the same modeling tools used for the 22 GHz H 2 
masers to model the 6.7-GHz methanol maser transition. 



4 RESULTS 

4.1 Maser distribution 

The resulting velocity integrated image of the 6.7-GHz methanol 
masers around Cepheus A HW2 is shown in Fig.Q] No maser emis- 
sion was detected at < 5 Jy beam -1 . The 12 maser features that were 
identified are shown in Fig.|4]and presented in TableQjwith position 
offsets (Aa and AS), velocity (Vlsr), peak flux (7 P eakX fractional lin- 
ear polarization (/>/) and polarization angle (x). 

The 6.7-GHz methanol masers make up a nearly perfect ellip- 
tical distribution around HW2 and its disc, and the maser distribu- 
tion observed with MERLIN matches that seen in European VLBI 
Network (EVN) and Japane se VLBI Network (JV N) observations 
(Torstensson et al, in prep.; Sug ivama et alj|2008l) . A least square 
fit of an ellipse to the MERLIN data indicates that the masers oc- 
cur in a ring with a radius of ~650 AU, a thickness of ~300 AU, a 
projection angle on the sky of = -78° and an inclination angle of 
i = 71°. The maser ring appears tilted with respect to the molecular 
(R = 580 AU, cf> = -56° and i = 68°) and dust disc (R = 330 AU, 
(p = -59° and i = 56°). Observations at high resolution (9x7 mas) 
with the EVN further indicate that the maser velocities are best 
fit by infall towards HW2 at ~1.7 km s' dTorstensson et"al]|2008l 
Torstensson et al., in pre p.). While the molecular d isc rotates with 
a velocity of ~5 km s' dJimenez-Serra et alj|2007l) . no significant 
rotation is observed in the maser region. 

4.2 Maser polarization and turbulent line-width 

We used the radiative transfer method described in § [3] to deter- 
mine the emerging brightness temperature, T b AQ, and intrinsic 
line-width, AV,-, of the observed maser features. As a direct fit re- 
quires a regular maser spectrum that does not suffer from blending, 
we were able to fit to 7 of the masers around Cepheus A HW2. 
The results of the fits are shown in Figs. [2] and [3] and given in 
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Figure 2. Results of the full radiative transfers-model fits, for those maser features from Fig. [4] that do not suffer from blending and for which linear 
polarization was detected. The fits yield the emerging maser brightness temperature Ti,ACl and intrinsic maser line-width AV,-. Contours indicate the significance 
intervals Ay 2 = 0.25, 0.5, 1, 2, 3, 7, with the thick solid contour indicating the l<x area. 



Table 1. Observation Results 



Feature 


Aa 


AS 


Vlsr 


^peak 


Pi 


X 




(arcsec)" 


(arcsec)" 


(km s 


) (Jy beam - 1 ) 


(percent) 


n 


1 


0.826 


0.010 


-4.14 


136.1 


0.15 ± 0.02 


37 ±6 


2 


-0.113 


0.038 


-2.78 


15.5 


< 0.75 




3 


-0.489 


0.161 


-1.82 


101.5 


0.40 ± 0.06 


-79 ±2 


4 


-0.528 


0.189 


-2.30 


131.1 


1.10 ±0.08 


-65 ± 10 


5 


-0.564 


0.212 


-2.60 


136.3 


1.01 ±0.03 


-56 ±2 


6 


-0.593 


0.228 


-2.73 


74.8 


1.88 ±0.05 


-56 ±4 


7 


-0.710 


0.270 


-3.68 


15.3 


1.48 ±0.07 


-80 ±2 


8 


-0.739 


0.308 


-3.65 


22.0 


1.02 ±0.07 


-61 ± 3 


9 


-0.805 


0.356 


-3.79 


32.1 


1.04 ±0.15 


-65 ± 1 


10 


-0.830 


0.430 


-3.70 


16.2 


0.41 ± 0.02 


-49 ±2 


11 


-0.364 


0.687 


-4.66 


11.2 


< 0.84 




12 


-0.379 


0.845 


-4.55 


8.8 


1.65 ±0.34 


-3 ±4 


8 Position offsets with respect to a 


= 22 ; '56 m 17.98 s , 6 = 62°01'49.390" 





Table [2] They indicate an emerging brightness temperature in the 
range 10 8 < TbAQ. < 10 9 K sr, implying that the masers are not 
saturated, as the ratio of stimulated emission rate over decay rate 
0.04 < R/T < 0.4, and that our assumption for AO. in § [3] is 
valid. It is worthwhile to note that our model calculations of the 
T b AQ assume a maser decay rate of T ~ 1 s~' . However, the model 
emerging brightness temperature scale linearly with T and conse- 
quently the saturation level (R/T) does not change. The intrinsic 



maser line-width lies in the range 1.0 < AV, < 1.5 km s _1 , with 
a weighted average AV; = 1.3 ± 0.2 km/s, which is much larger 
than the observed line-width of < 0.3 km s~' and shows that re- 
broadening, as a result of maser saturation, has not occurred yet. As 
the temperature ne eded to excite the 6.7-GHz masers is ~ 200 K 
dCragg et alj|2005h . the thermal line-width is ~ 0.7 km s The 
intrinsic line-width is thus significantly broader due to turbulence 
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Figure 3. The angle 6 between the maser propagation direction and the magnetic field vs. the fractional linear polarization P/. Panels are presented for the 
7 maser features for which an accurate fit could be made to the line profiles. The solid curves indicate the best fit emerging brightness temperatures and the 
dashed curves are the ±lcr intervals. The horizontal solid and dashed lines are the measured polarization fraction and the corresponding errors. The final panel 
indicates the remaining 5 maser features with curves for typical emerging brightness temperatures. In this panel, the dashed horizontal lines are the 3<x upper 
limits for the maser features without measured polarization, as these limits provide strong constraints on 8. 



in the maser region and we estimate the turbulent line-width to be 
AV « 1.1 km $~\ 

Using the observed P; and the fitted emerging brightness tem- 
perature we were then able to determine the angle between the mag- 
netic field and the maser propagation direction 9. The quoted errors 
on 9 are the most compact 68 percent probability interval, deter- 
mined by analysing the full probability distribution function. For 
those maser features where no radiative transfer fit was possible, 
including the features without detected linear polarization, we have 
used the range of brightness temperatures as suggested by the fitted 
maser features. The inclination angles 9 as well as the angle on the 
plane of the sky S for the maser features are given in Table[2] The 
errors in <p B are determined from the polarization angle x> with an 



added uncertainty of 5° as described in §|3] As seen in Tableland 
shown in Fig. [3] 9 is constrained to 9 > 9 CIlt at better than lcr level 
for the all but three of the maser features. Of those three, the best 
fit 9 for features 3 and 10 are still larger than # cr i,. For the majority 
of the maser features the magnetic field is thus likely perpendicular 
to the measured polarization vectors. Only for feature 1, the best 
fit 9 is less than 9 alt , meaning that the magnetic field is more likely 
parallel to x- 
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Figure 4. Total intensity (I; solid) and linear polarization (dashed) spectra for the 12 maser features detected above a level of 5 Jy beam in the MERLIN 
observations. The linear polarization spectra have been multiplied by a factor 20. 
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Table 2. Fit Results 



Feature AV; T b AQ 





</>B 


(kms 4 ) (togKsr) 


o 


(°) 



1 


1 


4 , +0.04 


8 - 08 -M 


47!" 


37 ±11 


3 


1 


30 ;8:8? 

Ju -0.11 


8.20^ 


«1 


11 ±7 


4 








73_f 


25 ± 15 


5 


1 


4fi +0.07 


8 59+ 035 

°- J -0.08 


7^ + 16 
' J -10 


34 ±7 


6 


1 


J -0.08 


8.90_-| 


77 + >° 


34 ±9 


7 

8 


1 


^c+0.11 
J -0.11 


8.78^;« 


74+15 

' -9 
7 9+l3 


10±7 
29 ± 8 


9 








72+14 


25 ±6 


10 


1 


OQ+0-04 
uy -0.10 




65^ 3 


41 ±7 


12 


1 


05+0.05 

UJ -0.14 


8.90.2 


72+8 
/z -10 


87 ±9 



5 DISCUSSION 

5.1 Methanol masers and the discs of Cepheus A HW2 



As noted in § 14.11 the 6.7-GHz methanol masers lie in an ellip- 
tical distribution around the molecular and dust discs of Cepheus 
A HW2. The masers are found at 650 AU, outside the molecular 
disc with a radius of 580 AU. The maser distribution also appears 
~ 20° tilted with respect to the dust and molecular discs and shows 
no sign of rotation. This indicates that the masers are not part of 
the molecular disc itself. The te mperature needed to maintain the 
masing conditions, 7* ~ 200 K dCragg et alj|2005h . is slightly less 
than the T = 250 ± 30 K found in th e outer regions of the molec- 
ular disc dJimenez-Serra et al]|2009l) . Also, t he hydrogen numb er 
density in the maser region, n Hl ~ 10 cm - 3 dCragg et al.i r2005) is 
approximately three orders of magnitude less than that in the disc 
dJimenez-Serra et alj|2009f) . Finally, the masers exhibit signs of in- 
fall and their velocity range covers only ~ 3 km s _I , which indicates 
that it is unlikely that they are part of the Cepheus A HW2 outflow. 
We thus suggest that the maser, while not associated with the disc 
directly, probe material that is being accreted onto the disc from the 
surrounding medium. 

5.2 The magnetic field strength 

Dynamic range limits in our MERLIN observations did not enable 
us to directly determine the magnetic field strength on the individ- 
ual maser features around Cepheus A HW2 and we obtained a 3<r 
limit of ~ 20 mG for the line-of-sight magnetic field. However, the 
magnetic field strength in the methanol maser region o f Cepheus A 
was m easured using the Effelsberg 100-m telescope ( IVlemmings! 
2008). Through Zeeman-splitting observations, the line of sight 
field strength was determined to be By = Bcos(6>) = 8.1 ± 0.2 mG 
dVlemmingsll2008h . which is mostly constant across the full maser 
spectrum. As this value is a flux averaged magnetic field deter- 
mined with a single dish telescope, we need to determine the 
amount of flux recovered and the number of maser features de- 
tected in our higher resolution MERLIN observations, to evaluate 
if the field measurement can be taken to be representative of the ob- 
served maser features. In Fig. [5] we show the comparison between 
the MERLIN and Effelsberg spectra. It is clear that the shape of the 
maser spectrum has changed significantly between the two epochs, 
though all individual maser features found with MERLIN can be 
identified in the Effelsberg spectrum. The chang e of the relative 
flux o f the individual maser features was noted in lSugivama et al.l 
(2008), who monitored Cepheus A with the Yamaguchi 32-m tele- 
scope. Comparing our spectra with their figure 1, shows that the 
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Figure 5. Total intensity spectra of the MERLIN observations (solid line) 
and the Effelsberg observations (dashed line) that were used to deter- 
mine the mag netic field strength in the 6.7-GHz methanol maser region 
I Vlemmings 2008) 



MERLIN spectrum is remarkably similar to what was observed at 
August 12 2007, with no more than 10 percent of the flux resolved 
out at high angular resolution. The MERLIN observations recover 
almost all of the single dish maser flux, and all features in the Ef- 
felsberg observations can be identified at high angular resolution. 
We can thus conclude that the Effelsberg magnetic field measure- 
ment can be taken as representative for true line-of-sight magnetic 

field in the maser region. 

As discussed in IVlemrningsl j2008h , there remains an esti- 
mated 25 percent uncertainty on the magnetic field strength due 
to the unknown accuracy of the Zeeman splitting coefficient of the 
6.7-GHz methanol maser transition. In addition to the measurement 
error of 0.2 mG on the magnetic field determination we thus add an 
additional 2 mG uncertainty to the magnetic field strength when we 
determine the effect of the magnetic field around Cepheus A HW2. 



5.3 Magnetic field morphology 

The full 3-dimensional magnetic field orientation was inferred us- 
ing maser radiative transfer models described in §[3] and is shown 
in Fig. [6] With the exception of one of the maser features that orig- 
inates from the mid-plane of the molecular disc behind the outflow, 
the inferred magnetic field in the plane of the sky is perpendicu- 
lar to the molecular and dust discs with a rms weighted average 
angle <p B = 26 ± 12°. Moreover, the error weighted average angle 
between the magnetic field and the line of sight 6 = 73 ± 5° is con- 
sistent with the inclination of the molecular disc and the outflow. 
The overall magnetic field orientation angle corresponds closely 
to the mag netic field direction observed in the encompassing dust 
envelope dCurran & Chrvsostomoull2007h . The discrepant polariza- 
tion angle found on the mid-plane maser feature is either due to the 
twisting of magnetic field lines towards the disc, or due to Faraday 
rotation caused by the ionized outflow. While the methanol masers 
at 6.7-GHz are not significantly affected by interstellar Faraday ro- 
tation, an electron density of ~200 citT 3 along a short 100 AU path 
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through the ionized outflow with a few mG magnetic field can cause 
the observed polarization angle of the masers behind the outflow to 
rotate by more than 10°. 



5.4 The role of the magnetic field 

Since we have obtained the angle between the magnetic field and 
the line of sight, we can directly determine the magnetic field 
strength around the protostellar disc. For fiy = 8. 1 ± 0.2 ± 2 mG and 
9 = 73±5°, the absolute field strength |B| = 23+' mG. From this, we 
can calculate a number of parameters to asses the role of the mag- 
netic field during the massive protostellar collapse. The Alfvenic 
mach number m a = a V3/ V A , where a = AVI V8 In 2 and AV is the 
turbulent line-width. The fits to the radiative transfer models indi- 
cate that the intrinsic maser line-width is much larger than the ther- 
mal line-width expected in the maser region with a typical kinetic 
temperature 7* ~ 200 K. Assuming the intrinsic line-width to be 
broadened due by turbulent motions, we find that AV ~ 1.1 km s~' . 
In the methanol maser region with a hydrogen number density 
«h 2 ~ 10 9 cnr 3 , the Alfven velocity V A = \B\I ^4np = 2.05 kms" 1 . 
The measured infall motions, at 1.7 km s" 1 , are thus slightly sub- 
Alfvenic and m a = 0.4. At 200 K the sonic mach number m s = 1.1, 
and consequently, the ratio between thermal and magnetic energy 
8 = 2(m a /m s ) 2 = 0.27. The magnetic field dominates the energies 
in the high density protostellar environment probed by the masers. 

We can also calculate the mass to magnetic flux ra- 
tio M/fl) compared to the critical value of this ratio A = 
(M/®)/( M/3>) crit . The critical value (M/<D) C rii = Cft/V^, ~ 
0.12/ V<5 l lMouschovias & Spitzerlfl976l : iTomisaka et al.ll 19881) is 
taken from numerical models for the collapse of a spherical molec- 
ular cloud into a highly flattened structure and defines the criti- 
cal mass that can be supported against collapse by the magnetic 
field. When A < 1, the magnetic field prevents collapse, but when 
A becomes larger than 1 gravity overwhelms the magnetic field. 
In terms of hydrogen column density, A can be determined from 
A = 7.6 x IQ- 24 N(H 2 )/\B\, with the magnetic field in mG. Taking 
the maser region to be ~300 AU thick, as determined by our ob- 
servations, the hydrogen column density along the magnetic field 
lines is N(H 2 ) « 5 x 10 24 cnr 2 . This yields A = 1.7!^, indicating 
that the region is slightly super-critical, a condition needed for the 
collapse to proceed. Similar values for A were found using large- 
scale Z eeman-splitting observations of an ensemble of molecular 
clouds (Crutcher 1999). However, as those observations lacked the 
3-dimensional magnetic field information, the results were uncer- 
tain by of order a factor two, making it impossible to distinguish 
between sub- and super-critical cloud cores. 

Finally, comparing the methanol maser magnetic field strength 
with other field strength measurements in the high-density region 
around HW2 (n Hl > 10 5 cm" 3 ), w e find that t he field scales approx- 
imately as B oc 5 (Fig. 3 from Vlemmings 2008). This relation 
could be the consequence of ambipolar diffusion, but it also nat- 
urally occurs for the collapse of a spherical cloud with frozen in 
field lines, when the preferred infall direction is along the magnetic 
field. Such a collapse forms a flattened disc-like structure similar 
as observed around HW2. We thus suggest that the magnetic field 
plays a crucial role in regulating the final stages of the formation of 
the massive protostar Cepheus A HW2, as is commonly expected 
during low-mass star-formation. 



5.5 Future perspectives 

5. 5. 1 Cepheus A as a prototype methanol maser source ? 

Due to their distance and complexity, the relation between the 
methanol masers and the protostellar environment for most other 
maser sources is still unknown. However, Zeeman-splitting obser- 
vations have shown typical magnetic field strengths along the line 
of sight of ~1 2 mG in over 70 p ercent of the observed methanol 
maser sources (Vlemmings 2008). Strong and dynamically impor- 
tant magnetic fields are thus present in most massive star-forming 
region traced by methanol masers. Furthermore, a recent survey has 
revealed that ~30 per cent of the methanol ma ser sources display 
an elliptical structure (Bar tkiewicz et alj|2009f) . This suggests that 
the methanol masers of a large fraction of star forming regions are 
located in an interface between a disc or torus and the infalling ma- 
terial, similar to the inferred origin of the masers around Cepheus 
A. Our picture of magnetic field regulated infall and accretion to- 
wards Cepheus A HW2 thus potentially describes a large number 
of massive star-forming regions. 



5.5.2 Further e-MERLIN observations 

With the e-MERLIN upgrade jGarrington et alj|2004) . it will soon 
be possible to map multiple maser transitions simultaneously with 
the radio-continuum. By observing in full polarization mode, this 
will allow for a direct comparison between the magnetic field struc- 
ture determined from the masers and the continuum morphology. 
Single track observations will allow for mapping of massive star 
forming regions at different evolutionary stages down to a con- 
tinuum sensitivity approaching 3 /jJy beam -1 . Simultaneously, the 
sensitivity in a narrow 6.7-GHz maser frequency band will then 
be sufficient to detect linear polarization for > 75 percent of the 
star formin g regions with metha nol masers when assuming the flux 
analysis of Ivan der Waltl d2005l) . When analyzed in a similar way 
to the Cepheus A HW2 observations presented here, such observa- 
tions have the potential to significantly further our understanding 
of the magnetic field around massive protostars. 



6 CONCLUSION 

We have demonstrated the power of maser polarization observa- 
tion, and particularly those using methanol masers, in deducing 
the full 3-dimensional strength and structure of the magnetic field 
around massive protostars. The detection of a coherent magnetic 
field direction in maser features that individually have a size of 
a few AU, suggests that the masers do not probe isolated pock- 
ets of a shock compressed magnetic field. In that case, the polar- 
ization direction would be determined by the direction of com- 
pression for each maser individually and observing coherent po- 
larization angles would be unlikely. A similar conclusion could 
be drawn from recent observations of W75N at high angular res- 
olution with the EVN that reveal a uniform magnetic field direc- 
tion in individual methanol maser features spread over ~ 2000 AU 
dSurcis et "al] |2009h . The detection of a large scale magnetic field 
is further confirmed by the good agreement of the field direction 
with that determined using dust polarization observations of the 
entire Cepheus A region. The 3-dimensional magnetic field geom- 
etry around Cepheus A HW2 supports the theories in which mag- 
netic fields regulate the infall and outflow close to massive pro- 
tostars in a similar way as during low-mass star-formation, even 



Magnetic field regulated protostellar infall 9 




Figure 6. The 3-dimensional magnetic field structure around the massive protostar Cepheus A HW2. The top left panel corresponds to a viewing angle close 
to the observations while the three further panels have different viewing angles. Spheres indicate the masers, with the black vectors indicating the true magnetic 
field direction. The maser features are color-coded according to their velocity. T he molecular disc IPatel et al ]|2005h is indicated by the dashed line and light 
grey structure while the dust disc I Jimenez-Serra et al. 2007; Patel et al. 2005) corresponds to the dark grey ellipse. The blue- and red-shifted lobes of the 
collimated radio outflow <CurieletaL 2006) are also shown. The red lines indicate the proposed magnetic field morphology. 
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if the high-mass star-forming regions are considerably more com- 
plex. The strong magnetic field that appears to be threading the 
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